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Analytical Study of Catalytic Reactors 
for Hydrazine Decomposition 
Quarterly Progress Report No. 6 
October 15, 19^7 - January 14, 1968 
Contract NAS 7-458 


SUMMARY 


The Research Laboratories of United Aircraft Corporation under Contract 
NAS 7-458 vith the National Aeronautics and Space Administration are performing 
an analytical study of catalytic reactors for hydrazine decomposition. This 
report summarizes work performed during the seventh quarterly contract period 
from October 15, 19^7 to January l4, 1968 . Work during this reporting period 
has included the debugging and running of the computer program representing the 
two-dimensional steady-state model of a distributed-feed catalyzed hydrazine 
decomposition reaction chamber. Calculations have been made of the effects on 
steady-state temperature and reactant concentration distributions of nonuniform 
radial injection and of catalyst bed configurations exhibiting both radial and 
axial nonuniformities. 

Empirical predictions have been developed of one-dimensional steady-state 
temperature and fractional ammonia dissociation profiles in hydrazine reactors 
packed with Shell 405 catalyst particles. The empirical correlations were 
developed on the basis of many runs made with the steady-state computer pro- 
gram developed during the first year of effort on the present contract. It 
was found that fractional ammonia dissociation and bulk fluid temperature are 
easily predicted for a broad range of operating conditions for cases in which 
most of the hydrazine decomposition occurs in the first few tenths of an inch 
of the reactor; this rapid hydrazine decomposition rate is associated with 
reactors packed with particles 25 mesh or smaller for approximately 0.2 in. 
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INTRODUCTION 


Under Contract NAS 7-458, the Research Laboratories of United Aircraft 
Corporation are performing analytical studies of the behavior of distributed- 
feed catalytic reactors for hydrazine decomposition. The specific objectives 
of this program are (a) to develop computer programs for predicting the tem- 
perature and concentration distributions in monopropellant hydrazine catalytic 
reactors in which hydrazine can be injected at arbitrary axial and radial 
locations in the reaction chamber and (b) to perform calculations using these 
computer programs to demonstrate the effects of various system parameters on 
the performance of the reactor. 

Progress previously reported in the first three quarterly reports 
(Refs. 1, 2 and 3) and in the first annual report (Ref. 4) included the develop- 
ment of computer programs which describe the steady-state and transient behavior 
of a reactor system in which complete radial mixing in the free-gas (or liquid) 
phase was assumed. Progress described in the fourth and fifth quarterly reports 
(Refs. 5 and 6) included an extension of the computer program describing the 
steady-state model to permit radial as well as axial variations in temperature 
and concentrations . The program contains a description of the turbulent dif- 
fusion of heat and mass in the interstitial phase along with heat and mass 
diffusion within the catalyst particles and between the particles and the inter- 
stitial phase. During this reporting period attention has been focused on 
debugging and running the computer program representing the two-dimensional 
steady-state model in order to evaluate steady-state reactor behavior. In 
addition, the one -dimensional steady-state program has been used to develop 
empirical correlations which predict axial temperature and fractional ammonia 
dissociation profiles in hydrazine reactors packed with Shell 405 catalyst 
particles . 
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DISCUSSION 


Effort during the seventh quarterly reporting period of Contract NAS 7-458 
has involved (a) debugging the computer program representing the two-dimensional 
steady-state model, (h) running this program to evaluate the effects on steady- 
state temperature and reactant concentration distributions of nonuniform radial 
injection and of catalyst bed configurations exhibiting both radial and axial 
nonuniformities, and (c) using the one-dimensional program to develop empirical 
correlations to predict axial temperature and fractional ammonia dissociation 
profiles in hydrazine reactors packed with Shell 405 catalyst particles . This 
effort is described in detail in succeeding sections of this report. 

Two-Dimensional Steady-State Program 


A series of calculations was made using the two-dimensional steady-state 
computer program in order to examine the effectiveness of the two-dimensional 
model and to evaluate the effects on system performance of nonuniform radial 
injection and of catalyst bed configurations exhibiting both radial and axial 
nonuniformities. The calculated results illustrated in Figs. 1 through 12 
refer to a reactor 3 in- in diameter into which liquid hydrazine is injected 
at the upstream end of the reactor only. For these calculations the upstream 
chamber pressure was taken as 100 psia and the hydrazine feed temperature was 
taken as 530 deg R. 

Axial temperature profiles at various radial locations are plotted in 
Fig. 1 for a case in which a radial nonuniformity in mass flow rate, G, is 
represented as a step function (see Fig. l). In this case the catalyst bed 
packing was taken to consist of 25-30 mesh catalyst particles for the first 
0.2 in. and l/8 in. x l/8 in. cylindrical pellets for the remainder of the bed. 
This configuration is referred to in the figures as the "standard bed con- 
figuration". Turbulent diffusion of heat, which tends to reduce radial tem- 
perature gradients, is more pronounced in the downstream end of the reactor. 

Here the catalyst particle size is larger, and both eddy conductivity and eddy 
diffusivity are directly proportional to particle size. The consequences of 
radial heat conduction are complicated somewhat by the simultaneous turbulent 
diffusion of mass. Higher temperatures are associated with more hydrazine 
decomposition; thus high temperature regions may lose heat by radial conduction, 
but may gain hydrazine from adjoining low temperature regions by radial dif- 
fusion of mass. Subsequent decomposition of this hydrazine may lead to even 
higher temperatures. For the case considered here, these combined effects lead 
to the temperature distribution shown in Fig. 1. For comparison purposes, the 
axial temperature profile corresponding to a radially uniform mass flow rate 
of 3-0 lb/ft^-sec is also plotted in Fig. 1. This is the average mass flow 
rate calculated by averaging the actual mass flow rate profile over the cross- 
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sectional area of the reactor. The mole-fraction distributions of hydrazine 
and ammonia associated with the temperature distribution shown in Fig. 1 are 
illustrated in Figs. 2 and 3 respectively. 

The results of similar calculations made for a continuously varying 
injection profile are shown in Figs. 4 through 6. The mass flow rate profile 
used for these calculations is plotted in Fig. 3- As in the first set of 
calculations, the average mass flow rate is 3-0 lb/ft^-sec. 

The effects on temperature and reactant concentration distributions of 
two catalyst bed configurations exhibiting both radial and axial nonuniformities 
are illustrated in Figs. 7 through 12. For both of these configurations the 
mass flow rate was taken as uniform at 3*0 lb/ft^-sec. The calculated temper- 
ature, mole fraction of hydrazine, and mole fraction of ammonia distributions 
corresponding to the bed configuration shown in Fig. 7 are plotted in Figs. 7 
through 9 respectively. Similar calculations corresponding to the bed con- 
figuration shown in Fig. 10 are plotted in Figs. 10 through 12. 

One-Dimensional Steady-State Program 


A series of runs were made with the one -dimensional steady-state computer 
program in order to develop empirical correlations to predict axial temperature 
and fractional ammonia dissociation profiles in hydrazine reactors packed with 
Shell 405 catalyst particles . Empirical correlations were developed on the 
basis of about 65 runs representing different combinations of mass flow rates, 
pressures and catalyst bed configurations. It was found that fractional ammonia 
dissociation and bulk fluid temperature can be predicted using the equations 

I - Fractional Ammonia Dissociation = <$> 


and 


Tj = 1020 + £ 0.075 CP/1000) ] j- + 1535 


where 


$ = (0.66) (G/Z)° 28 | (0.55 Q° 17 - 0.17 


(1000/P) 


0.22 


+ 0 


l7 } 


and z and a are expressed in ft, G in lb/ft^-sec, P in psia, and T-? in deg R. 
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These equations are illustrated in Figs. 13 and 14 respectively for cases in 
which most of the hydrazine decomposition occurs in the first few tenths of an 
inch of the reactor; this rapid hydrazine decomposition rate is associated 
with reactors packed with particles 25 mesh or smaller for approximately 0.2 
in. For these cases the correlations depicted in Figs. 13 and l4 work well 
for axial distances greater than one inch and for values of pressure, P, 

"between 10 and 1000 psia, mass flow rate, G, between 1.44 and l4.4 lb/ft^-sec 
(0.01 and 0.1 lb/ in. -sec respectively) and equivalent spherical radius, a, 
between 0.001 and 0.01 ft. For a reactor packed with small ( < 25 mesh) 
particles for the first few tenths of an inch and larger particles thereafter, 
the particle radius, a, refers to the larger particles. 

In Figs. 13 and l4, Rocket Research experimental data are plotted along 
with the empirical predictions and the results of sample cases run using the 
one-dimensional steady-state program. Values of fractional ammonia dissocia- 
tion obtained from the steady-state program are plotted for axial locations 
between 1 and 6 inches while values of bulk fluid temperature obtained from 
the program are plotted only for axial locations between 3 and 6 inches . 

It should be emphasized that these empirical correlations do not correctly 
predict the behavior of reactors in which hydrazine decomposition is slow, 
for example reactors which are uniformly packed with large catalyst particles, 
such as l/8 in. x l/8 in. cylinders. The correlations work quite well though for 
catalyst bed configurations consisting of 25-30 mesh particles for the first 
0.2 in. and l/8 in. x l/8 in. cylindrical pellets for the remainder of the bed. 
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LIST OF SYMBOLS 

a Radius of spherical particle, ft 

A p Total external surface of catalyst particle per unit volume of bed, ft" 1 

c j Reactant concentration in interstitial fluid, lb/ft3 

C p Reactant concentration in gas phase within the porous particle, lb/ft J 
C P * Equals c p - (c p ) s , lb/ft3 

C F Specific heat of fluid in the interstitial phase, Btu/lb - deg R 
C F Average specific heat of fluid in the interstitial phase, Btu/lb - deg R 

C s Specific heat of catalyst particle, Btu/lb - deg R 

Dj Diffusion coefficient of reactant gas in the interstitial fluid, ft^/sec 

D p Diffusion coefficient of reactant gas in the porous particle, ft 1 / sec 

F Rate of feed of hydrazine from distributed injectors into the system, 

lb/ft3-sec 

g c Conversion factor, (lbm/lbf) ft/sec 1 

G Mass flow rate, lb/ft 1 - sec 

h Enthalpy, Btu/lb 

A c Heat transfer coefficient, Btu/ft^-sec-deg R 

H Heat of reaction (negative for exotheiroic reaction), Btu/lb 

k c Mass transfer coefficient, ft/sec 

k Q Reaction rate constant, equals sec" 1 

K p Thermal conductivity of the porous catalyst particle, Btu/ft-sec-deg R 

M Molecular weight, lb/lb mole 

M Average molecular weight, lb/lb mole 
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n Order of decomposition reaction 

P Chamber pressure, psia 

Q Activation energy, Btu/lb mole 

r Radial distance, ft 

i' het Rate of (heterogeneous) chemical reaction on the catalyst surfaces, 
lb /ft 3 -sec 

r hom Rate of (homogeneous) chemical reaction in the interstitial phase, 
lb/ft3 -sec 

R Gas constant, equals 10,73 psia - ft3/lb mole - deg R, or, 

Radius of reactor 

t Time, sec 

t Actual time minus time required, under steady-state conditions, for 

liquid hydrazine to flow from the reactor inlet to the interface between 
the liquid-vapor and vapor regions, sec 

T Temperature, deg R 

Wj Weight fraction of reactant in interstitial phase 

x Radial distance from the center of the spherical catalyst particle, ft 

Z Axial distance, ft 

a Preexponential factor in rate equation 

ot p Intraparticle void fraction 

/3 Equals (Cp) s H Op / [k p (T p ) s 1 

Y Equals Q/R(T p ) s 

§ Interparticle void fraction 

€ Eddy diffusivity, ft^/sec 

X Eddy conductivity, Btu/ft-sec-deg R 
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H- Viscosity of interstitial fluid, lb/ft - sec 
Pj Density of interstitial fluid, lb/ft3 

p s Bulk density of catalyst particle, lb/f t3 
Subscripts 

F Refers to feed 

i Refers to interstitial phase 

j Axial distance index in transient model 
P Refers to gas within the porous catalyst particle 

S Refers to surface of catalyst particle 

Superscripts 

J Refers to chemical species 

L Refers to liquid at vaporization temperature 
V Refers to vapor at vaporization temperature 
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FIG. 2 


TWO - DIMENSIONAL STEADY - STATE DISTRIBUTION 
OF MOLE - FRACTION OF HYDRAZINE 
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FIG. 3 


TWO-DIMENSIONAL STEADY - STATE DISTRIBUTION 
OF MOLE - FRACTION OF AMMONIA 
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FIG. 4 


TWO-DIMENSIONAL STEADY -STATE TEMPERATURE DISTRIBUTION 



AXIAL DISTANCE, z - FT 





MOLE - FRACTION OF AMMONIA IN INTERSTITIAL PHASE 
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TWO - DIMENSIONAL STEADY - STATE DISTRIBUTION 
OF MOLE - FRACTION OF AMMONIA 
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TWO - DIMENSIONAL STEADY - STATE DISTRIBUTION 
OF MOLE - FRACTION OF HYDRAZINE 

P = 100 PSIA G = 3.0 LB/FT 2 - SEC 


CATALYST BED CONFIGURATION: 
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MOLE-FRACTION OF AMMONIA IN INTERSTITIAL PHASE 
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TWO - DIMENSIONAL STEADY - $ 
OF MOLE - FRACTION 

P = 100 PSIA G 

CATALYST BED CONF 


ALL 25-30 MESH 
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FIG. 9 
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FIG. 10 


TWO-DIMENSIONAL STEADY - STATE TEMPERATURE DISTRIBUTION 


P = 100 PSIA 


G = 3.0 LB/FT -SEC 


CATALYST BED CONFIGURATION: 



0.0167 FT 
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FIG. 12 


TWO - DIMENSIONAL STEADY - STATE DISTRIBUTION OF MOLE - FRACTION OF AMMONIA 

P= 100 PSIA G = 3.0 LB/FT 2 - SEC 

CATALYST BED CONFIGURATION: 



R = 0.125 FT 






FRACTIONAL AMMONIA DISSCIATION 
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21 FIG.14 

EMPIRICAL PREDICTION OF INTERSTITIAL 
GAS TEMPERATURE 


SYMBOL 

UNITS 

P 

PSIA 

Tj (i) 

DEG R 



O 0.2 0.4 0.6 0.8 1.0 

0 + [0.075 (P/1000)] 



